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Abstract 

In this article, we derive the scalar Klein-Gordon equation from the 
formal information theory framework. The least biased probability dis- 
tribution is obtained, and the scalar equation is recast in terms of a 
Fokker-Planck equation in terms of the imaginary time, or a Schroedinger 
equation for the proper time. This method yields the Green's function 
parametrized by the proper-time, and is related to the results of Schwinger 
and Dewitt. The derivation can then allow the use of potentials as con- 
straints along with the Hamiltonian or moments of the evolution. The 
information the- oretic, analogously the maximum entropy method, also 
allows one to ex- amine the possibility of utilizing generalized and non- 
extensive statistics in the derivation. This approach yields non-linear 
evolution in the Klein- Gordon-like partial differential equations. Fur- 
thermore, we examine the Klein-Gordon equation in curved space-time, 
and we compare our results to results obtained from path integral ap- 
proaches. 

1 Introduction 

In this article, we derive the scalar Klein-Gordon equation from within the 
information theory (or maximum entropy) framework. The least biased 
probability distribution is obtained from the maximization of the entropy 
with suitable con- straints. The probability then evolves under the scalar 
wave equation with a parametrized mass term, the Klein-Gordon equa- 
tion parametrized by the proper time. The scalar equation is recast as a 
Fokker-Planck equation in terms of the parametrization of imaginary time, 
or a Schroedinger equation for the proper time. This method is akin to 
the Feynman parametrization of the mass term to solve the Klein-Gordon 
wave equation and the derivation of the photon Green's function, as well 
as the Schwinger-Dewitt proper-time formalism. The derivation also al- 
lows one to use a higher dimensional version of the methods of quantum 
mechanics, specifically, the use of potentials in a linear superposition with 
the Hamiltonian of the evolution via the use of observable constraints in 
the entropy maximization procedure. We then use the method to examine 
the underlying stochastic differential equations and the microscopic evo- 
lution. We explore the solution to the wave equation for a particle moving 



in curved spacetime and given the moments obtained from the underlying 
stochastic microsocpic evolution, and obtain the short-time transitional 
(conditional) probability from the maximum entropy variational princi- 
ple. 

In the past, the proper-time formalism has gained acceptance in the 
treatment of the relativistic wave equation both in flat and curved space- 
time. The method simply involves a parametrization of the mass term, 
and recasts the Klein-Gordon equation into a partial differential equation 
(PDE) with an extra dimension (here) being that of the parametrization 
(the proper time). The solution is then obtained by transforming from the 
proper-time to the mass term [21 [3]. Also, well known analytic continua- 
tion methods allow us to recast the result in terms of a higher dimension 
Schroedinger equation with the subsequent connection to known methods 
of low energy quantum mechanics and the Schwinger-Dewitt proper time 
formalism. What has not been done is to consider the resultant partial 
differential equation as a diffusion-like equation in the context of probabil- 
ity theory and the maximum entropy method IS [2, or equivalently from 
within information theory |T . A mathematical reason for this is that the 
equations of the wave equation and Klein-Gordon form of equations are 
elliptical equations, and maximum entropy and information theory derive 
PDEs that are parabolic, of the diffusion, Fokker-Planck and Schroedinger 
forms, and that have 'friction'-like terms that allow for equilibriation with 
a heat bath in the thermodynamics perspective. 

This interpretation of the relativistic elliptical wave equation and Klein- 
Gordon equations as parametrized higher dimensional parabolic equations 
following the work of Feynman will allow us to recast the problem as one of 
maximum entropy and information theory. From derived least biased dis- 
tributions and their evolution parabolic PDEs we will obtain an equivalent 
description of the problem in terms of the underlying stochastic differential 
equations. For the curved-spacetime case this will mean two new formu- 
lations, an extensive statistics and the recently developed nonextensive 
statistics formulation of the stochastic gravity approach. The moments 
obtained from the stochastic description can then be used to obtain the 
short-time transition probability for the evolution. The Klein-Gordon 
partial differential equation is 

(1) 

and we rewrite the imaginary time as —it = r . We can obtain this 
equation from the maximum entropy method, or equivalently the infor- 
mation theoretic framework as the information measure is equivalent to 
the entropy with appropriate constants of proportionality relating the bits 
of information to the measure of disorder the entropy measured by the 
Boltzmann constant. We begin the derivation by stating that we know the 
moments of the variables {x, y, z, r) and maximize the Gibbs-Boltzmann 
entropy < S >= — J PlnP, equivalently minimizing the Shannon infor- 
mation measure, given the following constraints parametrized by A 
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Maximizing the entropy we obtain the least biased probability distri- 
bution function (PDF) normalized with N 



(6) 

This PDF solves the diffusion-like partial differential equation PDF 
dl'ix,!f, z, r.. A) _ Dr B'^P d'-P D„ D'-P D, ()'-P 

(7) 

We can rewrite the PDF via a Laplace or alternatively Fourier trans- 
form as 

and substitute in the PDE Eq. ([T| and obtain 



, a'- zj, a- D, a- d. a^ ,,,, 

(9) 

With the identification of the diffusion constants as Dx,y,z ~ 2 and 
Dr — 2c and (fi = — "^^^2 » we recover the Klein-Gordon equation as in 
Eq.([T}. We can also obtain a 5D Schroedinger equation using this method. 
The PDE is recast in this form if we take the time-like parameter and 
Wick rotate about the imaginary axis to A = —ia. The resulting PDF is 
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Schroedinger-like 



dI'{x..t/...-:,r,n) _ Dr a-I' D,.D'-I' D^O-r D.iS-P 



(10) 



and we can identify the diffusion constants as relating the 'mass' to the 
variance as Dx,y.z = = o"^ and Dr = — Alternatively, one 

can work with the Hamiltonian directly. From the Legendre transform 
of the relativistic Lagrangian, one can obtain the generalized canonical 
momenta, and write a relativistic Hamiltonian and any potentials as the 
observables to be used as constraints in the entropy maxmization proce- 
dure. Thus the constraints will be the < H >, and any potentials to 
be included such as external potentials or interactions < V >, and the 
entropy < 5* > is maximized to yield the least biased probability. In the 
case oi < V >— and flat space-time, the least biased probability will 
solve the massive Klein-Gordon equation as in Eq. 



2 Connection with nonextensive statis- 
tics 



Having obtained a maximum entropy derivation of the KG type equation, 
a question arises as to the form of the equations obtained if one were 
to use a generalized entropy. Instead of the Gibbs-Boltzmann entropy 
which yields the extensive statistics, one can use a generalized form of the 
entropy, the nonextensive entropy, which has become well known after 
the work of G. Tsallis O [7] and which yields non-linear partial differ- 
ential equations for the evolution of the distribution. In cosmology, the 
non-extensive statistics have been shown to model 'small' systems such as 
galactic matter distributions, where the range of the interaction (re: grav- 
itational) is on the same scale as the system size [TD], the solar neutrino 
problem [9] and systems that exhibit complex and non-linear dynamics 
. The present derived result should prove of interest as a new nonlinear 
theory as nonlinear Klein-Gordon models are of importance in high energy 
particle physics. The non-extensive statistical entropy, or analogously the 
incomplete information theoretic measure [l] is 



This form of entropy is known to yield power law distributions. Maxi- 
mizing the entropy with the observable moments as constraints yields the 
least biased power law distribution, 



1 - j r"{x,T,X)dxdr 



(11) 
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where Z^(A)/3(A) = const, and Z{X) is the partition function and is 
related to the normaUzation as usuaL This distribution solves the non- 
linear partial differential equation 



A) 



ax 



( h 



a^ 

dy^ ax^ 



]l'{£,r.X]- 



(13) 



where we have set h = c — 1 and absorbed any constants into the variables 
(a;, y, z, r). We can obtain the solution of the Klein-Gordon-like equation 
if we Laplace or alternatively Fourier transform according to 



(14) 



Some comments should be made regarding Eq.( 13 1. The nonlinear form of 



the PDE Eq.( 13 \ is one that by q-parametrization is a controllable devia- 
tion from the Gaussian solution of q— > 1 towards the power law solution 
of Eq. ( 12 I . The nonlinear terms of the traditional nonlinear Klein-Gordon 



equations of particle high energy physics are solved by functions that are 
deviations from the Gaussian, squeezed or stretched Gaussians. These 
nonlinear deviations are present by parametrized nonlinearity q and the 
PDE equation Eq.(13l should be considered as a new class of nonlin- 
ear Klein-Gordon equations. Another comment is to be made regarding 
curved space-time applications of this equation Eq.(13l. The curvature 
terms are included as the diffusion coefficients of the PDE Eq.Q as the 
curvature tensor terms are for diagonal tensors or tensors that can be 
transformed to diagonal tensors the inverses of the diffusion coefficients. 
These curvature terms cause the PDE Eq.Q to become nonlinear and 
the solutions of which become deviations from the Gaussian as discussed. 
Therefore the possibility exists here that the complicated nonlinearities 
due to nonlinear terms included to describe high energy particles and non- 
linear terms introduced to describe curved space-time can be described or 
modeled by the g- parametrized nonlinearity of Eq.(13l which has known 
solutions of the power-law form of Eq.( |12[ ). Also it is important to note 
that a recent advance in nonextensive quantum mechanics ^12^ is promis- 
ing for a relativistic Schroedinger-like PDE approach to nonextensive rel- 
ativistic quantum mechanics. This will be pursued in subsequent work. 



3 Stochastics 

The diffusion equation Eq.Q is a Fokker-Planck equation. This macro- 
scopic evolution equation for the PDF allows us to obtain the underlying 
Ito-Langevin stochastic differential equations [J. The individual trajec- 
tories are 

fix = ^U^iiivix),, (hj = /u;r/n ■ ( A) ^ 

<ii = sfDldW{X]., dT = </D^dlV{X)r-, 
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where the dWx,y,z,T are Wiener processes, and are composed of Gaus- 
sian white noise with delta correlations. The case of the non-extensive 
statistics derived Klein-Gordon-like PDE Eq.( |13[ ) is somewhat more com- 
plicated. The underlying stochastic evolution is of the statistical feedback 



form [6], and from Eq. ( 13 1 we obtain 



rf; = s/DJ'^-"{x.r,\)diV{X),, dr = ^/D , l'^-"{x,T, \)diV{X), . 



(16) 



We note again that the macroscopic distribution PDF appears in the un- 
derlying microscopic stochastic evolution. This statistical feedback 'loop' 
between the macroscopic evolution of probability distribution and the 
microscopic stochastic trajectory(s) can be viewed as the source of the 
nonlinearity in the PDE from this perspective. The modification of the 
trajectories via coupling to the large scale macroscopic evolution, alterna- 
tively the memory effect or the historical trends squeezes the PDF from the 
usual Gaussian type into a power-law, and can be a prototype for nonlin- 
ear deviations from the usual parametrized Klein-Gordon based Gaussian 
statistics. We note again that due to the presence of curvature, metrics 
with nonlinear terms will contribute terms to the Klein-Gordon equation 
that will cause deviations from a Gaussian [2] and that can be modeled by 
the q-parametrized nonlinearity. We will examine the regular form of this 
evolution in the context of curved spacetime in the next section. Mapping 
the regular curved space-time parametrized evolution to the nonextensive 
statistics g-parametrized evolution we leave for subsequent work. 



4 Klein-Gordon Equation in Curved Space- 
time 

The Klein-Gordon equation in curved spacetime from within the proper- 
time formalism of Schwinger and Dewitt has been treated, at least from 
the Gaussian approximation, by Bekenstein [21 13] . The equation is deriv- 
able from an action principle 



(17) 



The Euler-Lagrange equation then yields the Klein-Gordon equation in 
curved space-time, with c = Ti = 1 

{^"df^n,, - III" -iHH = a 

(18) 

We will solve the general case first, with R = R(x). The solution 
method will depend on the use of stochastic calculus and the maximum 
entropy method. We define the two-point function P{x,t,\;x ,t' ,\'). 
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This function evolves in Eq. (181. Writing generally, the backwards PDE 
is 

BP B'' 



3l = V'^' 



, ...V V- + h;—-(,B. /X^^T, A , r', A'). 
BX \ flr- J ^^g^ 

We can transform this equation into a standard form if we transform 
according to 



r) = / ;5(f, r) = / . ^ tfjy', 
The transformed equation is 

ai' ( 3'' a'' cP \ ,^ .d^i' 

dX \da^ BiP 3'!^; Bt^ 



(20) 



(21) 



We redefine the function P as P{d, r, A) = 7re ^'"^ . Upon substitution 
and gathering terms, we have 

Btt _ f B- B^ B- \ .-iB'ir r,\ '^^ 
, .BU Btt , ,BU Btt d^'lJ d'^lJ Br'U 

(22) 

The curvature term is eliminated if we make the identification of the 
first term with ^R. Upon replacing the original arguments, we have 



+ (-2 fc-^ — — + (-2 6:! — —. 

" Oil a» - at at , , 

■■^ (23) 

This equation is now in the backwards Chapman-Kolmogorov form. The 
two-point function, or the conditional distribution, is a solution of this 
equation, as it is a solution of the forward Chapman-Kolmogorov equation, 
the Fokker-Planck equation. 



-|-[^„{x)t]-|-[^;(.-?)^j . 

(24) 
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where wo have renamed the drift coefficients as Fx{x) = 2b'^x^, and 
similarly for the other coefHcients and have dropped the prime superscripts 
of the variables. The underlying stochastic differential equations are 



dx = h\ ix]dX + ^/bl{x\T)dl\\ (A), dr = ^/f^^^rjdW^ (A) , 

(25) 

where the dW are Wiener processes and are delta correlated as before. 
The time averages of the descretized stochastic processes are {Ax = x—Xo) 

(Ajt) =^-,(^')AA, (Ai^) =hHS', t')AA, 

(At) = 0, {AT=)=i?(f',T')AA. ^2^^ 

These averages coincide with conditional expectation values such as 

(Aj:|f',T', X) = J AaT5r(if, ir, iA|f , t", A')iiif<iAT = F4i')AX, 
(Ai'I^'.t'.A) = / ijT'sr(A£,AT,AAii',T',A')dAiiJAT = 4j(/,T')AA, 

^ (27) 

and so on for the other coordinate conditional expectation values. We 
solve for the short-time conditional (transition) probability distribution 
as follows. The conditional variances about the means are now the con- 
straints as we maximize the conditional entropy 

S{x',t', X') — - j jr(AE, At, AX^S', t' ,X')lnn(As, At, AAjS', r' ,X')dA£dAT. 

(28) 

Explicitly we maximize the following expression 

SliS)]+S [g, {(ir - {{&,, - {&y)f)] + - (A^))=) + {At=)] 

(29) 

The majcimization yields the least biased probability 

ir(Ax,Ar,AA|f ,t',A') = ^_e-(<'-('i'-(A'))'+f,(Ai'-(Ai())'+ti.(A— {As))=+;f,.i. =1 
Z{X] 

(30) 

The Lagrange multipliers are obtained from the partition function via 
functional relationsliipK l)ctween the multipliers and the moments. As an 
example, the x-coordinate multipliers are 

and the partition function is Z = J n. Upon integration, this yields 

/2F 12^ [2Tt phr 
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and we can solve for the multipliers using the relationships Eq.(31l to 
obtain 

1 1 



(33) 



and similarly for the multipliers in the other coordinates, in terms of 
their moments. Gathering the terms together and substituting for the 
moments, we have 



/'(Af, Ar, A r ,A ) = 

. t' IlIA ^. 21j^,| F' . t' iiA ^, 2 | F' . t ' 1 i A p 2 | F' , t ' I l1 A 



V'3?rft2.(J'. T']AX ^/27!f>l[x', T'lAA x/a^rfc^C:?, t') AA v'2rffcf {^'. t ') AA ' 

(34) 

We can compare this result with the proper-time formalism of Schwinger- 
Dewitt and the results of Bekenstein [21 13] if we wick rotate to the proper 
time A — —is. The Green's function from the path integral approach is 



6'(AjT Ar, A) = i{-lrrj«) ^ f-;'^^^ A5"[l + fits + /■)(<«)" + ...l 

■ ■ " ■ (35) 

The factors in the expression are (y{x; x) = o[x'; x) and stands for half the 
proper distance squared dS^ = g^,^dx^dx,^, or minus half the proper time 
squared, and contains the metric coefflcients, which we have re-written 
as inverses and are contained in the diffusion coefficients x,y,z,T ■ The 
quantity A is A{x;x') = —g{x)2 DvMg{x')2 , and Dvm is the VanVleck- 
Morreta determinant Dvm ~ Det{— g^^g^,^ ) • The flat space-time result is 

G(Aa;, At, A)o = i(47ris)~ ss^ A2, and Bekenstein solves the Gaus- 
sian case for the geodesic extreme paths and obtains 

(?(Ai',AT. Al = «(lT(,s)~'p''^p-'^'-^l''"A(jtf')^. 

(36) 

We note again that the diffusion coefficients are the inverse of the metric 
coefficients and are here included in the factor a. The normalizations are 
expectedly different, as is the R curvature scalar term, as the solution in 
the present case takes into account terms higher in the expansion. The 
short-time transition probability is accurate to order 6'(AA). We regain 
the flat space-time case where J? = if we let the means vanish (the drift 
coefficients), as the scalar curvature is included in the x, y, z, r coordinates' 
means, and assume constant (unit) diffusion coefficients. 
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5 Conclusion 



In this article wc liavc derived the scalar Klein-Gordon equation from 
within themaximumentropymethod,obtaining the (Wick rotated) proper- 
time analog and solutions. We have also generalized the approach to 
include other than Gibbs-Boltzmann entropies, namely the non-extensive 
statistics entropy of C. Tsallis. This resulted in a new nonlinear form of 
the Klein-Gordon equation, with its nonlinear proper time PDE. We then 
examine the Klein-Gordon equa- tion in curved space-time, and obtain a 
general solution in which the diffusion coefRents are related to the inverse 
metric coefHcionts. Also, the curvature scalar is then included in the 
drift. We compare our result with the results of Bekenstein and Dewitt. 
We leave the application to specific curved space-time metrics to future 
work. 
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